The SpRY Cas9 variant release the PAM sequence constraint for genome editing in the model plant Physcomitrium patens.

Genome editing via CRISPR/Cas has enabled targeted genetic modifications in various species, including plants. The requirement for specific protospacer-adjacent motifs (PAMs) near the target gene, as seen with Cas nucleases like SpCas9, limits its application. PAMless SpCas9 variants, designed with a relaxed PAM requirement, have widened targeting options. However, these so-call PAMless SpCas9 still show variation of editing efficiency depending on the PAM and their efficiency lags behind the native SpCas9. Here we assess the potential of a PAMless SpCas9 variant for genome editing in the model plant Physcomitrium patens. For this purpose, we developed a SpRYCas9i variant, where expression was optimized, and tested its editing efficiency using the APT as a reporter gene. We show that the near PAMless SpRYCas9i effectively recognizes specific PAMs in P. patens that are not or poorly recognized by the native SpCas9. Pattern of mutations found using the SpRYCas9i are similar to the ones found with the SpCas9 and we could not detect off-target activity for the sgRNAs tested in this study. These findings contribute to advancing versatile genome editing techniques in plants.

Regulation of developmental gatekeeping and cell fate transition by the calpain protease DEK1 in Physcomitrium patens.

Calpains are cysteine proteases that control cell fate transitions whose loss of function causes severe, pleiotropic phenotypes in eukaryotes. Although mainly considered as modulatory proteases, human calpain targets are directed to the N-end rule degradation pathway. Several such targets are transcription factors, hinting at a gene-regulatory role. Here, we analyze the gene-regulatory networks of the moss Physcomitrium patens and characterize the regulons that are misregulated in mutants of the calpain DEFECTIVE KERNEL1 (DEK1). Predicted cleavage patterns of the regulatory hierarchies in five DEK1-controlled subnetworks are consistent with a pleiotropic and regulatory role during cell fate transitions targeting multiple functions. Network structure suggests DEK1-gated sequential transitions between cell fates in 2D-to-3D development. Our method combines comprehensive phenotyping, transcriptomics and data science to dissect phenotypic traits, and our model explains the protease function as a switch gatekeeping cell fate transitions potentially also beyond plant development.

A DELAY OF GERMINATION 1 (DOG1)-like protein regulates spore germination in the moss Physcomitrium patens.

DELAY OF GERMINATION 1 is a key regulator of dormancy in flowering plants before seed germination. Bryophytes develop haploid spores with an analogous function to seeds. Here, we investigate whether DOG1 function during germination is conserved between bryophytes and flowering plants and analyse the underlying mechanism of DOG1 action in the moss Physcomitrium patens. Phylogenetic and in silico expression analyses were performed to identify and characterise DOG1 domain-containing genes in P. patens. Germination assays were performed to characterise a Ppdog1-like1 mutant, and replacement with AtDOG1 was carried out. Yeast two-hybrid assays were used to test the interaction of the PpDOG1-like protein with DELLA proteins from P. patens and A. thaliana. P. patens possesses nine DOG1 domain-containing genes. The DOG1-like protein PpDOG1-L1 (Pp3c3_9650) interacts with PpDELLAa and PpDELLAb and the A. thaliana DELLA protein AtRGA in yeast. Protein truncations revealed the DOG1 domain as necessary and sufficient for interaction with PpDELLA proteins. Spores of Ppdog1-l1 mutant germinate faster than wild type, but replacement with AtDOG1 reverses this effect. Our data demonstrate a role for the PpDOG1-LIKE1 protein in moss spore germination, possibly alongside PpDELLAs. This suggests a conserved DOG1 domain function in germination, albeit with differential adaptation of regulatory networks in seed and spore germination.

JGI Plant Gene Atlas: an updateable transcriptome resource to improve functional gene descriptions across the plant kingdom.

Gene functional descriptions offer a crucial line of evidence for candidate genes underlying trait variation. Conversely, plant responses to environmental cues represent important resources to decipher gene function and subsequently provide molecular targets for plant improvement through gene editing. However, biological roles of large proportions of genes across the plant phylogeny are poorly annotated. Here we describe the Joint Genome Institute (JGI) Plant Gene Atlas, an updateable data resource consisting of transcript abundance assays spanning 18 diverse species. To integrate across these diverse genotypes, we analyzed expression profiles, built gene clusters that exhibited tissue/condition specific expression, and tested for transcriptional response to environmental queues. We discovered extensive phylogenetically constrained and condition-specific expression profiles for genes without any previously documented functional annotation. Such conserved expression patterns and tightly co-expressed gene clusters let us assign expression derived additional biological information to 64 495 genes with otherwise unknown functions. The ever-expanding Gene Atlas resource is available at JGI Plant Gene Atlas (https://plantgeneatlas.jgi.doe.gov) and Phytozome (https://phytozome.jgi.doe.gov/), providing bulk access to data and user-specified queries of gene sets. Combined, these web interfaces let users access differentially expressed genes, track orthologs across the Gene Atlas plants, graphically represent co-expressed genes, and visualize gene ontology and pathway enrichments.

Improved prime editing allows for routine predictable gene editing in Physcomitrium patens.

Efficient and precise gene editing is the gold standard of any reverse genetic study. The recently developed prime editing approach, a modified CRISPR/Cas9 [clustered regularly interspaced palindromic repeats (CRISPR)/CRISPR-associated protein] editing method, has reached the precision goal but its editing rate can be improved. We present an improved methodology that allows for routine prime editing in the model plant Physcomitrium patens, whilst exploring potential new prime editing improvements. Using a standardized protoplast transfection procedure, multiple prime editing guide RNA (pegRNA) structural and prime editor variants were evaluated targeting the APT reporter gene through direct plant selection. Together, enhancements of expression of the prime editor, modifications of the 3' extension of the pegRNA, and the addition of synonymous mutation in the reverse transcriptase template sequence of the pegRNA dramatically improve the editing rate without affecting the quality of the edits. Furthermore, we show that prime editing is amenable to edit a gene of interest through indirect selection, as demonstrated by the generation of a Ppdek10 mutant. Additionally, we determine that a plant retrotransposon reverse transcriptase enables prime editing. Finally, we show for the first time the possibility of performing prime editing with two independently coded peptides.

Genome-wide specificity of plant genome editing by both CRISPR-Cas9 and TALEN.

CRISPR and TALENs are efficient systems for gene editing in many organisms including plants. In many cases the CRISPR-Cas or TALEN modules are expressed in the plant cell only transiently. Theoretically, transient expression of the editing modules should limit unexpected effects compared to stable transformation. However, very few studies have measured the off-target and unpredicted effects of editing strategies on the plant genome, and none of them have compared these two major editing systems. We conducted, in Physcomitrium patens, a comprehensive genome-wide investigation of off-target mutations using either a CRISPR-Cas9 or a TALEN strategy. We observed a similar number of differences for the two editing strategies compared to control non-transfected plants, with an average of 8.25 SNVs and 19.5 InDels for the CRISPR-edited plants, and an average of 17.5 SNVs and 32 InDels for the TALEN-edited plants. Interestingly, a comparable number of SNVs and InDels could be detected in the PEG-treated control plants. This shows that except for the on-target modifications, the gene editing tools used in this study did not show a significant off-target activity nor unpredicted effects on the genome, and did not lead to transgene integration. The PEG treatment, a well-established biotechnological method, in itself, was the main source of mutations found in the edited plants.

Physcomitrium patens Protoplasting and Protoplast Transfection.

Protoplast production with the moss Physcomitrium (Physcomitrella) patens has a long and successful history. As a tool, it has not only been the base of reverse genetic studies covering research fields as diverse as development, metabolism, or gene network regulation but also allowed its development as a bioengineering platform for protein production. We present here a standardized protocol for protoplast production from Physcomitrium (Physcomitrella) patens protonemata. Additionally, we detail procedures for their transfection, their plating for optimal regeneration, and three alternative selection approaches. To improve the consistency of protoplast regeneration, we describe a new option for protoplast embedding. The use of an alginate matrix to regenerate moss protoplast alleviates the use of warm agarized medium. Thus, it optimizes transformed protoplast survival without any morphological detrimental effect or impact on transfection efficiency.

Prime Editing in the model plant Physcomitrium patens and its potential in the tetraploid potato.

Since its discovery and first applications for genome editing in plants, the clustered regularly interspaced short palindromic repeats (CRISPR)-Cas9 technology has revolutionized plant research and precision crop breeding. Although the classical CRISPR-Cas9 system is a highly efficient tool for disruptive targeted mutagenesis, this system is mostly inefficient for the introduction of precise and predictable nucleotide substitutions. Recently, Prime Editing technology has been developed, allowing the simultaneous generation of nucleotide transitions and transversions but also short defined indels. In this study, we report on the successful use of Prime Editing in two plants of interest: the plant model Physcomitrium patens and the tetraploid and highly heterozygous potato (Solanum tuberosum). In both cases editing rates were lower than with other CRISPR-Cas9 based techniques, but we were able to successfully introduce nucleotide transversions into targeted genes, a unique feature of Prime Editing. Additionally, the analysis of potential off-target mutation sites in P. patens suggested very high targeting fidelity in this organism. The present work paves the way for the use Prime Editing in Physcomitrium patens and potato, however highlighting the limitations that need to be overcome for more efficient precision plant breeding.

A vertically transmitted amalgavirus is present in certain accessions of the bryophyte Physcomitrium patens.

In the last few years, next-generation sequencing techniques have started to be used to identify new viruses infecting plants. This has allowed to rapidly increase our knowledge on viruses other than those causing symptoms in economically important crops. Here we used this approach to identify a virus infecting Physcomitrium patens that has the typical structure of the double-stranded RNA endogenous viruses of the Amalgaviridae family, which we named Physcomitrium patens amalgavirus 1, or PHPAV1. PHPAV1 is present only in certain accessions of P. patens, where its RNA can be detected throughout the cell cycle of the plant. Our analysis demonstrates that PHPAV1 can be vertically transmitted through both paternal and maternal germlines, in crosses between accessions that contain the virus with accessions that do not contain it. This work suggests that PHPAV1 can replicate in genomic backgrounds different from those that actually contain the virus and opens the door for future studies on virus-host coevolution.

Single Nucleotide Polymorphism Charting of P. patens Reveals Accumulation of Somatic Mutations During in vitro Culture on the Scale of Natural Variation by Selfing.

Introduction: Physcomitrium patens (Hedw.) Mitten (previously known as Physcomitrella patens) was collected by H.L.K. Whitehouse in Gransden Wood (Huntingdonshire, United Kingdom) in 1962 and distributed across the globe starting in 1974. Hence, the Gransden accession has been cultured in vitro in laboratories for half a century. Today, there are more than 13 different pedigrees derived from the original accession. Additionally, accessions from other sites worldwide were collected during the last decades. Methods and Results: In this study, 250 high throughput RNA sequencing (RNA-seq) samples and 25 gDNA samples were used to detect single nucleotide polymorphisms (SNPs). Analyses were performed using five different P. patens accessions and 13 different Gransden pedigrees. SNPs were overlaid with metadata and known phenotypic variations. Unique SNPs defining Gransden pedigrees and accessions were identified and experimentally confirmed. They can be successfully employed for PCR-based identification. Conclusion: We show independent mutations in different Gransden laboratory pedigrees, demonstrating that somatic mutations occur and accumulate during in vitro culture. The frequency of such mutations is similar to those observed in naturally occurring populations. We present evidence that vegetative propagation leads to accumulation of deleterious mutations, and that sexual reproduction purges those. Unique SNP sets for five different P. patens accessions were isolated and can be used to determine individual accessions as well as Gransden pedigrees. Based on that, laboratory methods to easily determine P. patens accessions and Gransden pedigrees are presented.

Challenges of in vivo protein localization in plants seen through the DEK1 protein lens.

During the last 25 y, fluorescent protein tagging has become a tool of choice to investigate protein function in a cellular context. The information gathered with this approach is not only providing insights into protein subcellular localization but also allows contextualizing protein function in multicellular settings. Here we illustrate the power of this method by commenting on the recent successful localization of the large membrane DEK1 protein during three-dimensional body formation in the moss Physcomitrella patens. But as many approaches, protein tagging is not exempt of caveats. The multiple infructuous (failed) attempts to detect DEK1 using a fluorescent protein tag present a good overview of such potential problems. Here we discuss the insertion of different fluorescent proteins at different positions in the PpDEK1 protein and the resulting unintended range of mutant phenotypes. Albeit none of these mutants generated a detectable fluorescent signal they can still provide interesting biological information about DEK1 function.

Characterisation of evolutionarily conserved key players affecting eukaryotic flagellar motility and fertility using a moss model.

Defects in flagella/cilia are often associated with infertility and disease. Motile male gametes (sperm cells) are an ancestral eukaryotic trait that has been lost in several lineages like flowering plants. Here, we made use of a phenotypic male fertility difference between two moss (Physcomitrella patens) ecotypes to explore spermatozoid function. We compare genetic and epigenetic variation as well as expression profiles between the Gransden and Reute ecotype to identify a set of candidate genes associated with moss male infertility. We generated a loss-of-function mutant of a coiled-coil domain containing 39 (ccdc39) gene that is part of the flagellar hydin network. Defects in mammal and algal homologues of this gene coincide with a loss of fertility, demonstrating the evolutionary conservation of flagellar function related to male fertility across kingdoms. The Ppccdc39 mutant resembles the Gransden phenotype in terms of male fertility. Potentially, several somatic (epi-)mutations occurred during prolonged vegetative propagation of Gransden, causing regulatory differences of for example the homeodomain transcription factor BELL1. Probably these somatic changes are causative for the observed male fertility defect. We propose that moss spermatozoids might be employed as an easily accessible system to study male infertility of humans and animals in terms of flagellar structure and movement.

DEK1 displays a strong subcellular polarity during Physcomitrella patens 3D growth.

Defective Kernel 1 (DEK1) is genetically at the nexus of the 3D morphogenesis of land plants. We aimed to localize DEK1 in the moss Physcomitrella patens to decipher its function during this process. To detect DEK1 in vivo, we inserted the tdTomato fluorophore into PpDEK1 gene locus. Confocal microscopy coupled with the use of time-gating allowed the precise DEK1 subcellular localization during 3D morphogenesis. DEK1 localization displays a strong polarized signal, as it is restricted to the plasma membrane domain between recently divided cells during the early steps of 3D growth development as well as during the subsequent vegetative growth. The signal furthermore displays a clear developmental pattern because it is only detectable in recently divided and elongating cells. Additionally, DEK1 localization appears to be independent of its calpain domain proteolytic activity. The DEK1 polar subcellular distribution in 3D tissue developing cells defines a functional cellular framework to explain its role in this developmental phase. Also, the observation of DEK1 during spermatogenesis suggests another biological function for this protein in plants. Finally the DEK1-tagged strain generated here provides a biological platform upon which further investigations into 3D developmental processes can be performed.

PEATmoss (Physcomitrella Expression Atlas Tool): a unified gene expression atlas for the model plant Physcomitrella patens.

Physcomitrella patens is a bryophyte model plant that is often used to study plant evolution and development. Its resources are of great importance for comparative genomics and evo-devo approaches. However, expression data from Physcomitrella patens were so far generated using different gene annotation versions and three different platforms: CombiMatrix and NimbleGen expression microarrays and RNA sequencing. The currently available P. patens expression data are distributed across three tools with different visualization methods to access the data. Here, we introduce an interactive expression atlas, Physcomitrella Expression Atlas Tool (PEATmoss), that unifies publicly available expression data for P. patens and provides multiple visualization methods to query the data in a single web-based tool. Moreover, PEATmoss includes 35 expression experiments not previously available in any other expression atlas. To facilitate gene expression queries across different gene annotation versions, and to access P. patens annotations and related resources, a lookup database and web tool linked to PEATmoss was implemented. PEATmoss can be accessed at https://peatmoss.online.uni-marburg.de.

Analysis of the Localization of Fluorescent PpROP1 and PpROP-GEF4 Fusion Proteins in Moss Protonemata Based on Genomic "Knock-In" and Estradiol-Titratable Expression.

Tip growth of pollen tubes, root hairs, and apical cells of moss protonemata is controlled by ROP (Rho of plants) GTPases, which were shown to accumulate at the apical plasma membrane of these cells. However, most ROP localization patterns reported in the literature are based on fluorescent protein tagging and need to be interpreted with caution, as ROP fusion proteins were generally overexpressed at undefined levels, in many cases without assessing effects on tip growth. ROP-GEFs, important regulators of ROP activity, were also described to accumulate at the apical plasma membrane during tip growth. However, to date only the localization of fluorescent ROP-GEF fusion proteins strongly overexpressed using highly active promoters have been investigated. Here, the intracellular distributions of fluorescent PpROP1 and PpROP-GEF4 fusion proteins expressed at essentially endogenous levels in apical cells of Physcomitrella patens "knock-in" protonemata were analyzed. Whereas PpROP-GEF4 was found to associate with a small apical plasma membrane domain, PpROP1 expression was below the detection limit. Estradiol-titratable expression of a fluorescent PpROP1 fusion protein at the lowest detectable level, at which plant development was only marginally affected, was therefore employed to show that PpROP1 also accumulates at the apical plasma membrane, although within a substantially larger domain. Interestingly, RNA-Seq data indicated that the majority of all genes active in protonemata are expressed at lower levels than PpROP1, suggesting that estradiol-titratable expression may represent an important alternative to "knock-in" based analysis of the intracellular distribution of fluorescent fusion proteins in protonemal cells.

The Physcomitrella patens gene atlas project: large-scale RNA-seq based expression data.

High-throughput RNA sequencing (RNA-seq) has recently become the method of choice to define and analyze transcriptomes. For the model moss Physcomitrella patens, although this method has been used to help analyze specific perturbations, no overall reference dataset has yet been established. In the framework of the Gene Atlas project, the Joint Genome Institute selected P. patens as a flagship genome, opening the way to generate the first comprehensive transcriptome dataset for this moss. The first round of sequencing described here is composed of 99 independent libraries spanning 34 different developmental stages and conditions. Upon dataset quality control and processing through read mapping, 28 509 of the 34 361 v3.3 gene models (83%) were detected to be expressed across the samples. Differentially expressed genes (DEGs) were calculated across the dataset to permit perturbation comparisons between conditions. The analysis of the three most distinct and abundant P. patens growth stages - protonema, gametophore and sporophyte - allowed us to define both general transcriptional patterns and stage-specific transcripts. As an example of variation of physico-chemical growth conditions, we detail here the impact of ammonium supplementation under standard growth conditions on the protonemal transcriptome. Finally, the cooperative nature of this project allowed us to analyze inter-laboratory variation, as 13 different laboratories around the world provided samples. We compare differences in the replication of experiments in a single laboratory and between different laboratories.

The Physcomitrella patens chromosome-scale assembly reveals moss genome structure and evolution.

The draft genome of the moss model, Physcomitrella patens, comprised approximately 2000 unordered scaffolds. In order to enable analyses of genome structure and evolution we generated a chromosome-scale genome assembly using genetic linkage as well as (end) sequencing of long DNA fragments. We find that 57% of the genome comprises transposable elements (TEs), some of which may be actively transposing during the life cycle. Unlike in flowering plant genomes, gene- and TE-rich regions show an overall even distribution along the chromosomes. However, the chromosomes are mono-centric with peaks of a class of Copia elements potentially coinciding with centromeres. Gene body methylation is evident in 5.7% of the protein-coding genes, typically coinciding with low GC and low expression. Some giant virus insertions are transcriptionally active and might protect gametes from viral infection via siRNA mediated silencing. Structure-based detection methods show that the genome evolved via two rounds of whole genome duplications (WGDs), apparently common in mosses but not in liverworts and hornworts. Several hundred genes are present in colinear regions conserved since the last common ancestor of plants. These syntenic regions are enriched for functions related to plant-specific cell growth and tissue organization. The P. patens genome lacks the TE-rich pericentromeric and gene-rich distal regions typical for most flowering plant genomes. More non-seed plant genomes are needed to unravel how plant genomes evolve, and to understand whether the P. patens genome structure is typical for mosses or bryophytes.

Combination of the Endogenous lhcsr1 Promoter and Codon Usage Optimization Boosts Protein Expression in the Moss Physcomitrella patens.

The moss Physcomitrella patens is used both as an evo-devo model and biotechnological production system for metabolites and pharmaceuticals. Strong in vivo expression of genes of interest is important for production of recombinant proteins, e.g., selectable markers, fluorescent proteins, or enzymes. In this regard, the choice of the promoter sequence as well as codon usage optimization are two important inside factors to consider in order to obtain optimum protein accumulation level. To reliably quantify fluorescence, we transfected protoplasts with promoter:GFP fusion constructs and measured fluorescence intensity of living protoplasts in a plate reader system. We used the red fluorescent protein mCherry under 2x 35S promoter control as second reporter to normalize for different transfection efficiencies. We derived a novel endogenous promoter and compared deletion variants with exogenous promoters. We used different codon-adapted green fluorescent protein (GFP) genes to evaluate the influence of promoter choice and codon optimization on protein accumulation in P. patens, and show that the promoter of the gene of P. patens chlorophyll a/b binding protein lhcsr1 drives expression of GFP in protoplasts significantly (more than twofold) better than the commonly used 2x 35S promoter or the rice actin1 promoter. We identified a shortened 677 bp version of the lhcsr1 promoter that retains full activity in protoplasts. The codon optimized GFP yields significantly (more than twofold) stronger fluorescence signals and thus demonstrates that adjusting codon usage in P. patens can increase expression strength. In combination, new promotor and codon optimized GFP conveyed sixfold increased fluorescence signal.

An intragenic mutagenesis strategy in Physcomitrella patens to preserve intron splicing.

Gene targeting is a powerful reverse genetics technique for site-specific genome modification. Intrinsic homologous recombination in the moss Physcomitrella patens permits highly effective gene targeting, a characteristic that makes this organism a valuable model for functional genetics. Functional characterization of domains located within a multi-domain protein depends on the ability to generate mutants harboring genetic modifications at internal gene positions while maintaining the reading-frames of the flanking exons. In this study, we designed and evaluated different gene targeting constructs for targeted gene manipulation of sequences corresponding to internal domains of the DEFECTIVE KERNEL1 protein in Physcomitrella patens. Our results show that gene targeting-associated mutagenesis of introns can have adverse effects on splicing, corrupting the normal reading frame of the transcript. We show that successful genetic modification of internal sequences of multi-exon genes depends on gene-targeting strategies which insert the selection marker cassette into the 5' end of the intron and preserve the nucleotide sequence of the targeted intron.